A renewed look at rf in medium 



(N 

o 

(N 

u 

Or 

< 



O 

=5 



> 

O 

l> 

en 
O 
(N 



% 



Youngshin Kwon, 1,2 >u Su Houng Lcc, 2 '|j] Kcnji Morita, 3 '[*| and Gyorgy Wolf 4 '!! 

1 Department of Physics, Soongsil University, Seoul 156-743, Korea 

2 Department of Physics and Institute of Physcis & Applied Physics, Yonsei University, Seoul 120-749, Korea 

3 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

4 Wigner RCP, RMKI, H-1525 Budapest, POB 49, Hungary 

We revisit the question of whether the Ua(1) symmetry is effectively restored in hot and dense 
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I. INTRODUCTION 

The breaking of the Ua(1) symmetry is an operator 
relation that remains valid even when the spontaneously 
broken chiral symmetry is restored. However, whether 
its effect on the rf mass survives even when chiral sym- 
metry is restored is a phcnomenological question that has 
caught the interest of many researchers [H-Q ■ The ques- 
tion has recently been revived as the RHIC data on two 
pion Bose- Einstein correlation at y/s = 200 GeV Au+AU 
collision seems to suggest the quenching of the rf mass in 
medium [7H9|. Its partial quenching in nuclear medium 
is also of great interest as such effects could be probed in 
finer detail in nuclear target experiments [10h12| . 

The restoration of Ua(1) symmetry will depend on two 
important ingredients; its relation to chiral symmetry 
breaking and the effects from topological configurations. 
How the former and latter contribute to the n-point cor- 
relation functions have been clarified in references jj] and 
[a, Q. Combining both effects, it was shownQ that in 
the chiral limit, with Nf flavors, the symmetry will effec- 
tively be restored in correlation functions composed of up 
to Nj — 1 points [l3|. This means, for example, that for 
Nf = 3, Ua(1) symmetry will be effectively restored in 
the two-point functions when chiral symmetry is restored. 
Still, the argument is based on correlation functions and 
does not explicitly relate the rf mass to the other pseudo- 
scalar masses. To establish this relation, we revisit the 
Witten-Veneziano (WV) formula [IJ, LL5[ for the rf mass 
in vacuum and generalize it to finite temperature. Al- 
though the formula is obtained in the large N c limit, we 
will obtain an explicit relation that relates the rf mass to 
condensates and two-point correlation functions at finite 
temperature. Therefore the symmetry restoration pat- 
tern observed in the two-point function will be reflected 
in the rf mass. For the physical case of Nf = 3 with ex- 
plicit quark masses, this result implies that the mass of 
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rf will become degenerate with the other pseudo-scalar 
mesons up to the quark masses when chiral symmetry is 
restored. 

The paper is organized as follows. In section[Hj we will 
revisit the main results in ref. [5j. We will then review 
and generalize the WV formula in section Hill Wc discuss 
the rf mass in Sec. IV. The summary will be given in the 
last section. 



II. CORRELATION FUNCTIONS AND THE 

u A (i) symmetry 

Here, we start with a brief summary of the main re- 
sult given in ref.[5|. The starting point is the Euclidean 
partition function of QCD: 



Z[J] 



D[A]er SYM T)ct[IJ) + m q ] 



E z t J ]- 



(1) 



where Sym = (1/4)F 2 . The second line writes the parti- 
tion function in terms of topological configurations with 
the topological charge v = (g 2 /32) J d A xFF. The whole 
integral in the first line is a positive definite quantity 
[J], which we will denote as d\i for later convenience. 
The topological configurations are always accompanied 
by n+(nJ) number of right-handed (left-handed) fermion 
zero modes such that v = n+ — n_. In such topologi- 
cally non-trivial configurations, the fermion determinant 
comes with special chirality such that partition function 
with v = I can be written more explicitly as follows: 



Z[J\u=x = f D[A]u=ie- s ^ 



r Det'[^- 



xdet 



d A x4> (x)m q il) (x) 



(2) 



where the prime in the fermion determinant means that 
the chiral zero modes rpo have been explicitly taken out 
into the second determinant. Therefore, in the chiral 
limit m q = 0, the topological configuration does not con- 
tribute to the partition function as the fermion deter- 
minant gives zero. However, these terms do contribute 



in the correlations functions and select out the rj' from 
the other pseudo-scalars. Higher topological configura- 
tions will contribute at higher point functions when there 
are sufficiently many external legs to saturate the zero 
modes. 

To see this, consider a two point function of a generic 
quark bilinear: 

n T (x)=(q(x)Tq(x),q(0)Tq(0)) 

1 



d/i 



Tr[S A (x,0)TS A (p,x)T] 



Ti[Sa(x,x)T]Ti[Sa(0,0)T] + (zero mode) 



(3) 



the first, second and third term being the connected, dis- 
connected terms and possible zero mode contribution re- 
spectively. Sa is the quark propagator in the presence 
of the gauge field. When V contains a flavor matrix, the 
contributions from the disconnected diagrams are identi- 
cally zero. 

The results in refs.|4[ and Q can be summarized as 
follows. When one takes the difference between the two- 
point functions of chiral partners, the difference vanishes 
when chiral symmetry is restored. When the difference 
is taken between those composed of currents that are 
related by a chiral transformation and an extra Ua(1) 
transformation, there will be an extra contribution from 
the zero modes. As an example, the difference between 
a pseudo-scalar and rf is given as follows in SU(2): 



U 7r (x)-U n '(x) = — l d\i 



Ty[S a (x,x) 15 }Ti[Sa(0,0) 15 ] 



7? / d Vv= 



4Mx)^o (x)M0)M0) 



(4) 



In ref. [4| , it was shown that the first term goes to zero in 
the chiral limit when chiral symmetry is restored. This 
result is in fact independent of the number of flavors and 
also valid when the 75 inside the trace is replaced by 
other gamma matrices such as 7^ or 7^75. 

The zero-mode contributions appearing in the second 
term of Eq.Q come from the topological configuration 
in Eq. ([2]) and are responsible for the appearance of the 
11^(1) effect. However, when Nf > 2, the second de- 
terminant in Eq.© will have 2Nf zero mode lines and 
hence the zero- mode contributions in Eq.Q will be pro- 
portional to 0(m q t ) and vanish in the chiral limit. 
Therefore, when chiral symmetry is restored, the Ua(1) 
breaking effect will not appear in the two-point func- 
tions. However, this does not necessarily mean that the 
rf mass will become degenerate with the other pseudo 
scalar mesons because the coupling of the currents to the 
rf might just go to zero. Therefore, let us look at a re- 
lation that directly relates the mass of 77' to the chiral 
order parameters. 



III. WITTEN-VENEZIANO FORMULA 
A. WV formula at zero temperature 

As a first step of this study, we review the derivation 
of the WV mass formula [14[ . We start with the gluonic 
correlation function defined in a pure glue theory: 



U(k) 



d 4 xe l 



'(TGG(x)GG(O)). 



(5) 



One should note that, in the large iV c limit {16j|. Eq. ([5]) 
scales as order N% . There is also a well known low energy 
theorem for the correlation function at zero external mo- 
mentum U(k = 0) ^ 0, whose value we will come back 
in the next section. 

However, when massless quarks are added to the the- 
ory, the low energy theorem leads to the vanishing cor- 
relation function Ui q {k = 0) = 0, where the subscript Iq 
means the presence of light quarks, through the anomaly 
relation that relates the pseudo-scalar gluon current to 
axial current £ g <9 m <77m75<7 = N f j^GG, where G° v = 
l/2e /i!/a/ gG^o. This seems a little odd, because quark ef- 
fects are suppressed in large N c , but for the low energy 
theorem, the leading N c effect seems to be canceled by a 
sublcading iV c effect. The answer to this question led to 
the WV formula. 

In terms of the physical states, the correlation func- 
tions looks as follows when light quarks are added. 



u k (k) = -j2 



|(0|GG|ra th glucball)|' 
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In the spectral form, Eq. ((SJ) , the first term in the right 
hand side indicates contributions from glueballs, while 
the second term shows those from the meson composed 
of light quarks. One can show that the residue of the 
first term is of order N% whereas the quark effects are of 
order iV c [lB: 



(0 1 GG I n th glueball) = N c a r , 

(0\GG\ 



,th 



meson) = \/N r c„ 



(7) 



Since all the meson masses should have a smooth large N c 
limit O(l), the terms of quark effects, U\, are suppressed 
in 1/N C as expected. The resolution to the seemingly 
inconsistent result is by noting the existence of r/'-meson 
whose mass scales as order 1/N C . One then recovers the 
consistent low energy theorem Ui q (0) — Uq(0) + Ui(0) = 
0, if the second term is saturated by the ry'-meson and 
scales as N%. From this condition, one finds, 



17(0) = t/ (0) = - 



l(0|GG|r/)| 2 



N c 



(8) 



By using the V(1)a anomaly, 



{0\GGW) = --^(Q\d,j£W) 
a Nf 

4tt 1 



(9) 



a N f 



Nfm v ,U, 



Eq. ((5J becomes as follows: 
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(10) 



where Nf is the number of light flavors. In Eq. (|9]), we 
made use of /,,/ = f n to lowest order in N c . Eq. (fTUI) is 
the celebrated WV formula. 



B. WV formula at finite temperature 

Consider the correction to Eq. (fTUf at finite tempera- 
ture. As mentioned before, the correlation function in 
Eq. (|5|) is order N 2 , as can be seen by the two loops rep- 
resenting two gluon lines in Fig.[l]-(a). At finite tempera- 




fa) 





(c) 





(e) 



FIG. 1: Two sets of diagrams 

ture, the thermal correction could come from the thermal 
gluon or quark interactions. Fig.[T]-(b,c) show the ther- 
mal corrections to Uo(k) while Fig.[T]-(d,e) show those to 
U^k). 

The dominant thermal gluonic contribution to Uo(k) 
comes from Fig.[T]-(b) and scales as N 2 as in the vacuum 
scaling. The scaling comes as follows, 



where the factors are the same as before except for the 
last factor, which comes from the number of external 
quark lines, as can be seen from Fig.[T]-(c). Therefore, in 
the large N c limit, the thermal gluonic effect scales the 
same as in the leading vacuum scaling and will contribute 
to modifying Uo(k). 

As for the modification in the quark loops Ui(k), the 
thermal gluonic effects are shown in Fig.Q]-(d), and the 
thermal quark effects in Fig.[T]-(c). Both scale as 0(N C ), 
and can thus be neglected in the leading order correction. 

If the system is in the confined state, the hadronic side 
will be saturated by color singlet glucball, meson and nu- 
cleons. Here, the dominant contribution comes from the 
nucleons. One can show that the contribution from the 
nucleon to all figures in Fig.Q] scale as 0(N C ) because 
the nucleon contains N c quarks. On the other hand, the 
contributions from meson or glueballs are suppressed in 
l/Nc as the number of constituents arc finite. Hence, 
hadron effects can be neglected until near the phase tran- 
sition point where the density of states increases, after 
which one can use the quark and gluon degrees of free- 
dom. 

Therefore, same arguments hold as in the vacuum. 
Namely, the addition of quarks somehow has to cancel 
the leading N c behavior at k — 0. This cancelation can 
not be done by collective states, as quark collective states 
are also suppressed in large N c limit, and hence has to 
come from a modified rf contribution. All in all, a sim- 
ilar equation to Eq. ([SJ will hold at finite temperature, 
with |(0|GG|t/)| now defined at finite temperature at rf 
momentum zero. Moreover, it should be noted that the 
rf mass we are discussing now could be different from 
that of the pole mass as we are discussing the scalar part 
of the mass, which survives at k^ — > 0. A simplified 
example would be to assume that the small energy and 
momentum self energy has the following form, with a(T) 
and b(T) being the small corrections, 



E,/ = a{T)k 2 + b{T)k 2 + m 2 (T). 



The pole mass at 



y/m 2 + m 2 (T)/ y/1 - a(T). 



But 



would 
the mass we 



(13) 

be 

are 



talking about is \Jm 2 



\T). 



Nevertheless, Uq has a nontrivial correction at finite 
temperature as we will see in the following sections. 



(N c )x 



(N c ) 2 = Q(N 2 ). 



(11) 



IV. n' MASS AT FINITE TEMPERATURE 



which comes from the internal loop, two coupling and 
the number of external thermal gluons respectively. The 
contributions where the gluons couple directly to the cur- 
rents scale the same as in Fig.[T]-(b). On the other hand, 
the contributions from thermal quarks to Uo(k) scales as 
follows, 



A. Low energy theorem 

Now, Uo(0) can be obtained from the low energy the- 
orem. Here we use the derivation using the heavy quark 
expansion |18j . For technical reasons, we start from a 
slightly different definition of the correlation function. 



(A c )x 



1 



(A^) = 0(A C ), 



(12) 



P(k) = i 



d A xe l 



T 



^GG(x),^GG(0) 



(14) 



It can be shown [18[ that 

2 



P(k = 0) 



32tt 2 d(-l/4 ff g) \tt 



-G 



(15) 



Now, any matrix element with canonical dimension d 
should be proportional to the d th power of the scale 



Sir 2 
A = M exp - — 2 

b9o 



(16) 



with 6=11 — |JV/. Hence, the gluon condensate at finite 
temperature and density should be of the following form. 



- G2 > = A ^(x<x 

7T /T,/i \A A 



(17) 



where d = 4 and / is a generic function specifying the 
temperature and density dependence of the gluon con- 
densate. Then Eq. (flip]) becomes, 



P(k = 0) 



■!(<- r £-"£)<;°VM 



Now, combining Eq. ([5J and Eq. (fl~8]) . one finds, 

3 Q V |(0|GG|V)| 2 _2^ ^9 ^/V 

4n) m 2 ,, ~b\ dT ti dn)\n /t* 

(19) 
But now since we can make the identification of the left 
and right hand side only in the large N c limit, the right 
hand side should be calculated in the quenched approx- 
imation. This means that one should just read off the 
temperature dependence of the gluon condensate from 
the lattice calculation for pure gauge theory, and also 
take 6 = 11. Thus, the if mass is given by 



3a 



\(0\GG\ V ')\ 



'\|2 



(*-T&)<5&) 



(20) 



T,pure gauge 



B. Gluonic part 

In order to evaluate the in-medium if mass from 
Eq. (|2"0"j) , all we need are the temperature-dependence of 
the gluon condensate and the coupling of GG to if . First 
let us consider the denominator of Eq. (|2T)|) . It has been 
known for a long time, that the gluon condensate has 
contribution from the pcrturbative and non-perturbative 
contribution. Moreover, it was also known that at the 
critical temperature, the non-perturbative contribution 
changes abruptly, but does not vanish completely, and 
retains more than half of its non-perturbative value [19j- 

El. 

The effect of subtracting out the second term in the 
denominator of Eq. (f2TJ)) is to get rid of the perturbativc 
correction, or the seemingly scale breaking effect that is 
not related to scale breaking but due to the introduction 
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FIG. 2: T-dependence of gluon condensate and its derivative 
fitted to Wuppertal-Budapest lattice data in full QCD. 



of an external scale parameter T. The leading perturba- 
tive correction to the gluon condensate is proportional to 
g 4 (T)T 4 [Hi!. Therefore, assuming that the temper- 
ature dependence is of the following form, 



(-G 2 ) = G (T) + ag 4 T*, 

\7T IT 



(21) 



we find. 







i - T ^)(i; c;1 )r-{ d - T w) G ^- < 22 > 

if the temperature dependence of g is neglected |27|. The 
only temperature dependence that survives is Gq(T), 
whose scale dependence is coming from dimensional 
transmutation and not from the external temperature 
only. It is the non-perturbative part that dominates the 
behavior of the right hand side of Eq . (|20[) . 

In Fig.[5J we show the model fit 24| to the Wuppertal- 
Budapcst's full QCD data 25( together with the modified 
T-dependent gluon operator as it appears in Eq. (|2"0j) . 
In the quenched calculation, the only difference is that 
the change is takin g pl ace more abruptly near the new 
critical temperature [26| . One can see that the change in 
the denominator of Eq. ([2TJf will tend to reduce the if 
mass near the critical temperature. 



C. coupling to if 

The final step in obtaining the mass of if using 
Eq. (|20[) . when chiral symmetry is restored, is estimating 
the change of the coupling (0\GG\i)'). For that purpose, 
let us consider Ui q (k) in Eq. <j5j) in the full theory, but 
rewrite it in terms of the quark axial current using the 
anomaly relation. 

U(k) =i fd 4 xe lkx (TGG(x)GG(0)) 

2 r 

(Tgi7 M 7 5 g(x) <?i7„7 5 <7(0)) 



--k^k v i I d 4 xe ifc ' a 



-in 



aN 



(T<n M g Or) q'yM )) 



(23) 



where we have subtracted out the contribution from the 
conserved vector current. Using the previous terminol- 
ogy, when chiral symmetry is restored, the connected 
piece will cancel, as they are the same as the difference 



J 



between flavored chiral partners, and only the discon- 
nected pieces will remain. Assuming that the spectral 
sum starts from the 77', we find Eq. (|23|) can then be writ- 
ten as follows, 



U(k) 



\(0\GGW) 



k 2 



-^k^k v 



mi 



n 
4tt 



aN 



Tr[S , j4 (x,a;)i7 AI 75]Tr[S , A(0,0)z7^7 5 ] - Tr[S A (x,x)j l _ l }TT[SA(0,Oh v } 



(24) 



r 



However, the disconnected pieces are all of the same order 
in m q when chiral symmetry is restored. 

Tr[S A (x,x)] ~ Tr[5 A (x,x)r] ~ 0(m q ), (25) 

where T is a Hermitian gamma matrixQ. Since Eq. ([24]) 
is valid for any k, we find that 



(0\GGW) ~ OK), 



(26) 



when chiral symmetry is restored. Therefore, going back 
to Eq. ([20]) and making use of the previous discussions, 
we find that when chiral symmetry is restored, 



2 <«9>->-0 

m n , — > 0, 



(27) 



in the chiral limit. One concludes that in the large N c 
limit of QCD, rf mass will become degenerate with the 
other goldstone bosons. 



CONCLUSIONS 



mesons do not change their mass towards the phase tran- 
sition point, it is this extra U^(l) mass of rf that is going 
to be quenched in the chiral symmetry restored phase. 

Few remarks are in order. First, in the quenched ap- 
proximation, the changes of order parameters take place 
only near the phase transition point. This suggests that 
the effect of quenching might only be visible when the 
hadronization temperature is close to the phase transi- 
tion point as in the case of RHIC or LHC energies for 
example. Second, it is hard to make a quantitative es- 
timate on how much of this mass is quenched in the 
nuclear medium, as the effects of density are sublead- 
ing in the large N c limit. However, assuming Eq. ([20|) 
is an exact relation, we can use Eq. (f24| to approxi- 
mate (0\GG\rj') oc Tt[Sa(x,x)] oc (qq) and then use it in 
Eq. (|20|) to deduce m v i oc (qq) , assuming that the change 
in the gluon condensate is small in nuclear medium. 
Therefore, if the chiral order parameter reduces by 20% 
in nuclear medium the U^(l ) breaking part of the rf mass 
will also reduce by the same fraction. 



It should be noted that the rf mass that is being 
quenched is the part of the mass that comes from the 
breaking of the U^(l) symmetry. Going back to Eq. (fpQ|) 
and substituting the vacuum value of Eq. (|T5|) one finds, 



V 66 jx ir 



464 MeV, 



(28) 



where we have used U = 1 30 MeV and (fG 2 ) = 
(0.35GeV) 4 . This is smaller than the vacuum value of 
the r/' mass as expected. Assuming that the pseudo scalar 
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